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ABSTRACT 

The flow across a parallel sealing dam of the type that appears in 
noncontacting shaft face seals is analyzed for steady, laminar, subsonic, 
isothermal compressible flow with relative rotation of the sealing dam 
surfaces. Rotational effects on leakage rates, pressure profiles, and 
other physical quantities of interest are found. Conditions are given 
under which a hydrostatic analysis of the radial leakage flow is a valid 
approximation. The error in estimating the mass leakage by using the 
hydrostatic radial flow formula is given. Engineering formulas pre- 
sented include mass leakage flow rate, pressure and velocity distribu- 
tions, opening force and center of pressure. The physics of the flow 
are discussed. 

NOMENCLATURE 

A cross-sectional area, in.^ 

a speed of sound, ft/sec 

Cp specific heat at constant pressure, Btu/ (Ibm) (°R) 

C^ specific heat at constant volume, Btu/ (Ibm) (°R) 

D/Dt material derivative, 3/9t + u(3/3r) + (v/r)(3/30) +w(3/3z) 

E^( ) exponential integral function 
F sealing dam force, Ibf 


body force vector 
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film thickness 5 nominal, in. 

2 2 

dimensionless parameter = -3R2^ /SSffi 
sealing dam width, in. 


^ - 1 ) 


n=l 
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.n+1 


nn; 


Mach number 

mass flow, Ibm/min 

change in mass flow 

• 3 

dimensionless mass flow, 12pM£ T/irh P 2 
molecular weight of gas, Ibm/lb-mole 
an integer 

static pressure, psi 
pressure differential, psi 

smaller pressure of two pressure boundary conditions, psi 

2 

dimensionless pressure, P/PqU 
2 2 2 

dimensionless P , P /P 2 

universal gas constant, 1545.4 f t^lbf / (Ib^mole) (^R) ; 


R mean radius, (R^ + R2)/2, in. 

AR sealing dam length, R 2 - R^^, in. 

(it gas constant, R/m, f t-lbf / (Ibm) (^R) 

Re^ pressure flow Reynolds number in radial direction, pUh/p 

Re^ shear flow Reynolds number in circumferential direction, pRQh/\x 

* 

Re modified Reynolds number, Re^^h/AR 

r radial direction coordinate 
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r dimensionless radial coordinate, r/AR 

T temperature, 

U pressure flow reference velocity, ft/sec 

u velocity in r~direction or x-direction, ft/sec 

* 

u dimensionless velocity, u/U 

V reference shear flow velocity, ft/sec 

V velocity in 0-direction, ft/sec 

V dimensionless velocity, v/V 

w velocity in y-direction, ft/sec 

w dimensionless velocity, 

W reference velocity across film thickness, U(h/AR), ft/sec 

X transformed coordinate, t ^/^2 

center of pressure in radial or X-direction, in. 

X coordinate in pressure gradient direction 

y coordinate across film thickness 

y dimensionless coordinate, y/h 

z shear flow coordinate in Cartesian system 

V specific-heat ratio, G /C 

p V 

^ circumferential coordinate 

A second viscosity coefficient or coefficient of bulk viscosity 

2 

U absolute or dynamic viscosity, (Ibf ) (sec) /in. 

p 

V kinematic viscosity, ft'^/sec 

P density, (Ibf ) (sec^) /in.^ 

k 

P dimensionless density, P/Pq 

angular rotational velocity, rad/sec 

V Del operator, (3/8r)i + (l/r)(3/36)j + (3/3z)k 
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Subscripts t 


av average 

h based on film thickness 

r based on radius 

0 reference 

1 inner radius 

2 outer radius 

INTRODUCTION 

Some powerplantSj such as advanced jet engines, exceed the operat- 
ing limits of face contact seals (Refs* 1 and 2)* As a result, noncon- 
tact face seals are becoming necessary* A noncontact face seal which is 
pressure (force) balanced is shown in Fig* 1* In this seal the leakage 
and pressure drop occur across a narrowly spaced sealing dam, and the 
axial force associated with this pressure drop is balanced by a prede- 
termined hydrostatic closing force and a spring force* However, this 
configuration, with the sealing dam formed by parallel surfaces has an 
inherent problem; the force due to the pressure drop across the sealing 
dam is independent of film thickness; hence, there is no way of main- 
taining a preselected film thickness which will allow tolerable leakage 
and still have noncontact operation* Since the force is independent of 
film thickness, the design also lacks axial film stiffness for sufficient 
dynamic tracking of the stationary nosepiece with the rotating seal seat* 
The seal nosepiece must follow the seal seat surface under different 
operating conditions without surface contact or excessive increase in 
film thickness, which would yield high leakage* Some of these operat- 
ing conditions are axial runout, misalinement , and thermal deformation 
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(coning and dishing) , 

A successful method of maintaining a preselected film thickness 
and achieving axial film stiffness is to add self-acting lift pads, such 
as shrouded Rayleigh step bearings, to the conventional face seal (1 and 
2), This is illustrated in Fig. 2. The axial sealing dam force associ- 
ated with the pressure drop across the sealing dam, and the lift pad 
force, are balanced by the hydrostatic and spring closing forces. The 
gas bearing has a desirable characteristic whereby the force Increases 
with decreasing film thickness. If the seal is perturbed in such a way 
as to decrease the gap, the additional force generated by the lift pad 
will open the gap to the original equilibrium position. In a similar 
manner, if the gap becomes larger, the lift pad force decreases, and 
the closing force will cause the seal gap to return to the equilibrium 
position. Since a proper balance of the opening and closing forces must 
be found in order to determine a gap with a tolerable mass leakage, 
physical quantities of interest, such as pressure distribution and mass 
leakage, must be evaluated. 

In this paper only the sealing dam portion of the seal will be 
analyzed. The classical viscous, isothermal, subsonic, compressible 
flow analysis for this problem is well known (e.g., see Gross, (3)). 

The pressure distribution and mass leakage have been calculated for the 
parallel film hydrostatic case. Reference (4) analyzed the parallel 
film hydrostatic case including the effects of fluid inertia, viscous 
friction and entrance losses. Subsonic and choked flow conditions can 
be predicted and analyzed for both laminar and turbulent flows. Results 
showed good agreement with experiment. The effects of rotation of the 
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much greater than h^ the entrance region effects are neglected; hence^ 
the convective inertia forces are neglected. This means that the seal 
is treated as operating entirely in the viscous region » This case is 
contrasted to the case where the gap size is large ^ as illustrated in 
Fig, 4(b) for a radial diffuser. 6) The fluid film is isothermal. This 
means that all heat generated in the film is conducted away through the 
walls. This is a standard assumption of lubrication theory. The va- 
lidity of this assumption breaks down for cases of large thermal gra- 
dients in the sealing dam and when the frictional heating is high (e.g., 
small gap or high speed). However , thermal analysis of a gas film seal 
(2) shows that the sealing dam can be closely approximated by a constant 
temperature. The limit of the validity of isothermal duct flow analyses 
is that M < l//y (6). 7) The entrance Mach number is close to zero. 

This means that there cannot be a large axial flow source on one surface 
impinging on the radial surface, as in a hydrostatic bearing (Fig. 4(c)), 
8) The fluid velocity in the reservoir is considered to be negligible 
(stagnant), and thus its effects are neglected in this analysis. This 
model is, therefore, not valid for the radial step seal shown in 
Fig. 4(d) when radial velocities are significant. 

Governing Equations 

The analysis will be outlined briefly. Further details can be 
found in (7) , The governing flow equations for a compressible fluid 
with constant viscosity in vector notation are (8) 

Conservation of massi 

2^ + p? • V = 0 [1] 
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seal seat with respect to the seal nosepiece (sealing dam surfaces) on 
the radial pressure flow^ however, have been neglected in these com- 
pressible flow analyses* Snapp (5) has analyzed the effects of rotation 
for the incompressible flow case* 

This analysis was conducted 1) to investigate theoretically the 
effect of the relative rotation of the sealing faces on mass leakage, 
pressure and velocity distributions, net pressure opening force, and 
center of pressure (radial direction) for an isothermal, compressible, 
viscous flow with parallel sealing dam surfaces; 2) to determine under 
what conditions, if any, the hydrostatic radial flow analysis is suffi- 
cient for engineering accuracy* 

ANALYSIS 

Basic Model 

The sealing dam model, as shown in Figs. 3 and 4(a), consists of 
two parallel, coaxial, circular rings in relative rotation at a constant 
speed separated by a very narrow gap. A pressure differential exists 
between the rings’ inner and outer radii. The fluid velocities are 
small in both the inner -diameter cavity and outer -diameter cavity which 
bound the sealing dam. 

The model formulation is based on the following physical conditions 

1) The fluid is homogeneous, compressible, viscous, and Newtonian. 

2) The flow is steady and laminar (continuum flow regime) , and the body 
forces are negligible* 3) The bulk modulus is ignored (X = -2/3 u) . 

This is Stokes idealization (3) . This condition will be valid unless 
the gas is under high pressure, very dense (e,g., shock-wave structure), 
or rarefied* 4) The fluid behaves as a perfect gas. 5) Since AR is 
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Conservation of momentum (Navler-Stokes equations) 


DV 

Dt 


-VP - uV X (V X V) + ^|-Ui)v(V • V) + F 


[ 2 ] 


Equation of state for an isothermal process: 

P = P(p) [3] 

By using cylindrical coordinates (Fig. 3) and by applying the con- 
ditions assumed in the mathematical model, including the narrow gap 
approximation, the governing flow of a compressible fluid equations in 
dimensionless form reduce to 


Conservation of mass 
1 3 * * * 

— — (P r u ) 
r 3r 


0 


Conservation of momentum 
1) Radial direction 
* 2 * 


* 3u 

P U — 5 

3r 


-(i) 1 


* 2 * 
3P 1 3 u 

. * * *2 
Sr Re 3y 


2) Circumferential direction 


* 3v , * 3v , 

U — +W — + 

3r 3y 


* 

u V 


3) Axial direction 
* 

= 0, since 


3P 


3y 

Equation of state 
* 

P .... 

— r = constant 


3P 

3y 


2 * 

1 3 V 

* *2 
Re 3y 


(y 


[4] 


[5] 


[ 6 ] 


[7] 


[ 8 ] 


Equations [4] to [8] have been nondimens ionalized with the following 
variables; 
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Ik 
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r 


AR 


* 

y 




* 


w 


w 

W 




foUH 

y 


* 


u 




All the dimensionless terms with the exception of r are of unit mag- 

“k 

nitude (r is of greater than unit magnitude) . By neglecting the 

k 

entrance region and assuming Re << 1, the convective inertia forces 
are neglected* Hence the dimensionless radial momentum Eq. [5] be- 
comes 




2 * 
* V 

P ~ 

r 


9P 

9r 


2 * 

_ + -L- (Ml) 9 U 
* \h ) ^ * 


9y 


*2 


[9] 


It is important to note that if V/U < l^y^Re^(h/AR) , the radial 

pressure flow can be treated as uncoupled from the rotational shear flow* 

This condition justifies the use of radial hydrostatic leakage flow 

equations for gas film seals whose sealing faces are in relative rota- 

k 

tion. Again for Re « 1 Eq* [6] becomes 


2 * 
a V 


ay 


^2 


= 0 


[ 10 ] 


Now in dimensional form the equations to be solved are 
Conservation of mass (from Eq, [4])^ 


1^ (pru) = 0 [11] 

(This form of the continuity equation is not used but is replaced by the 
integrated form, which is shown to be the conservation of mass flow in 
the radial direction.) 
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Radial momentum (from Eq, [5])* 


PV^ _ I n 

r 9r . 2 

3y 


[ 12 ] 


(The centripetal inertia term on the left side is ignored in the compres- 
sible Reynolds lubrication equation,) 

Circumferential momentum (from Eq, [6])s 



[13] 


Now the above set of equations will be solved. Equation [13] may be in- 
tegrated directly. Applying the boundary conditions 

V = 0 at y = 0 V = rQ at y = h 

yields 


V 


h 


[14] 


Now Eq. [12] is solved by substituting for the circumferential ve- 
locity equation (Eq. [14]) and the equation of state (Eq. [8]). Since 
the circumferential and axial variations are shown to be negligible. 


P = P(r), and Eq. [12] becomes 

2 

a^u 1 dP ^0^ 2„ 

7T = 7d7-772-’^y 

3y P^h y 


[15] 


Integrating Eq. [15] in the y-direction and applying boundary conditions 
u“0 aty = 0 u = 0 aty = h 

yields 


u 


2ii dr 


- hy) - r 

12yPQh^ 



[16] 
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The mass flow at any radius is 


M = 2irrp 


u dy = - 


2 3 2 

,3 h pr P 

irh pr dP 0 

6p dr 20pP^ 


Substituting the perfect gas law (Eq. [8]) into Eq. [17] yields 


■ dp2 , 2^2 


.2 A •> 


dP / 3 


dr \ 5 P 


rP^ = - 


0_ li 
3 /r 


To facilitate the solution, transform Eq. [19] in the following manner: 


• * ^ -12pM £T 

3 2 

.rh^P^ 


\ = -TW- 


Thus, Eq. [19] becomes 


~ dP ~2~2 '* 

X^+K^XP =M 


The conservation of mass must be satisfied 


(which replaces the continuity Eq, [4]) or 


^ lx + K ^ (X^p^) = 0 
dX \ dX / dX 
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Equation [22] can be solved subject to the following boundary conditions; 





Equation [22] is a first-order nonhomogeneous ordinary differential 
equation with variable coefficients and can be solved by using the 
classical integrating factor method of solving differential equations. 
Applying the boundary conditions yields the following mass flow equation 
in dimensional form; 



[23] 


The integral in Eq. [23] can be expressed in terms of exponential inte- 
gral function (9). In this form Eq. [23] becomes 



The resulting squared pressure distribution has integrals which can also 
be written in terms of exponential integral functions. Hence 
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The total force per unit width is found from 


^ P^2 p 

L mxn 

Also the center of pressure is found from 


[26] 


X = 
c 


f^2-h 

vJq 


(P - dX 

min 

_ 

L 


[27] 


where is the smaller pressure of the two pressure boundary condi- 

tions. The integrations in Eqs. [23] to [27] are performed numerically 
in (10) . Additional equations of interest in designing gas film seals 
are also given in (10) . 

Equations [25] and [26] show that the opening force does not depend 
on film thickness. Thus there is no axial film stiffness for the par- 
allel sealing dam surface case^ and some auxiliary device ^ such as self- 
acting lift padSj must be used to obtain axial film stiffness (see 
Fig. 2). 

Circumferential Shear Flow 

Under the conditions of this analysis, which include many shaft 
face seal applications, the circumferential flow is the well known 
simple Couette shear flow. The circumferential velocity distribution 
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is then given by Eq. [14] which is 


^ h 


[ 14 ] 


The total power can be estimated by considering the viscous shear 
due to rotation only 


—2 2 
uR Q A 

Power = 

h 


[28] 


Armlytical Expressions for Special Cases 

Small Parameters 

For most practical seal problems j the value of the parameter -K^/2 
is less than 1, In this case the exponential integrals of Eqs • [24] 
and [25] can be evaluated by integrating an infinite series expansion. 
Thus analytical expressions can be found for the pressure distribution 
and mass leakage. Now, 



where 


[29] 



Thus, the mass flow Eq. [23] becomes 


[30] 
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The resulting squared pressure distribution equation becomes 



No Rotation 


For no rotation, = 0; hence, Z (r/R 2 ) =0, and after substitut- 
ing the perfect gas relation Eq. [31] becomes 


12y flT In 

h 

and Eq. [32] becomes 



[33] 


[34] 


Equations [33] and [34] are reducible to the form found in (3). In many 
gas film sealing dams AR is much less than R^. For this case 
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Eq, [33] becomes 


• - 4) 

12y ftTCR^ - 

and Eq<, [34] becomes 



[ 35 ] 


[36] 


This is reducible to the form of the narrow slot leakage equation; thus, 
hydrostatic radial flow can be approximated with a narrow slot (plane 
flow) analysis with small error if (R 2 ““ R^^) « 

For no rotation Eq, [27] can be readily integrated which results in 


1 

L 


2Pi(R2 


3 



- P . (R. - RJ 
min 2 1 


[37] 


The center of pressure can also be readily found from Eq. [28] which 
yields 



Error in Estimating Mass Leakage by Hydrostatic Radial-Flow Formula 

A formula will now be found to give the error that occurs in using 
the hydrostatic radial flow mass leakage equation when there is relative 
rotation of the sealing dam surfaces. Defining 
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Eq. [33] 


static 

and simplifying by applying conditions present in most sealing damsj 
namely ^ that 

1) For |k^/2| « 1, 


K, 


£ 1 + 


and 


exp 



1 + 


2R^ 


2) For » R 2 - Rj^, 


In 


AR 


R, 


3) £(R^/R2> = Kj^/4(R^/R2 - 1) 


The above expression yields; 

„2 


AR 


AM 


1 - 


2 2 
SR^n 

T£t“ 


1 - 



2 2 
P2 


M 


static 


R, 1 - 


4R . ^ 
R^ ^ 5KT 


- 1 


- 1 


[39] 


DISCUSSIOK 

Rotational Effects on Seal Leakage 

For laminar flow across an isothermal compressible flow sealing dam, 
the change in mass leakage due to relative surface rotation would appear 
to be important when the surface speed is very high. Such nn example. 


is a proposed space power generator (11) which requires small pressure 



differential j low leakage seals at high rotational speeds (36 000 rpm) • 
Figure 5 shows the effect of rotation on radial mass leakage for a pro- 
posed space power compressor seal* Both internally and externally 
pressurized modes are compared with the hydrostatic flow analysis. The 
film thickness varies from 0.1 mil to 0.3 mil which are practical oper- 
ating gaps. The sealed fluid is a helium-xenon mixture. The operat- 
ing conditions are? sealed pressure , 32 psia; ambient pressure , 

31.5 psia; sealed gas temperature, 330° F; inner radius, 3,90 in.; 
outer radius, 4.00 in. The radii used in this example are larger 
than in (11), From Fig, 5, the hydrostatic analysis predicts about a 
42% higher leakage than the rotational analysis when the seal operates 
in the externally pressurized mode. On the other hand, for the in- 
ternally pressurized mode the rotational analysis shows that leakage is 
about 27% higher than predicted by the hydrostatic analysis. Rotational 
effects would be expected to be important for this case since the param- 
eter V/U > l/C* (V/U = 310, = 10). 

The approximate formula given by Eq. [39] can be used to determine 
whether rotation is important in radial mass leakage calculations. If 
rotational effects are shown not to be important, the simpler hydro- 
static flow formulas can be used. Figure 6 illustrates that there is 
excellent agreement between the approximate formula (Eq* [39]), which 
includes only the first term of the convergent series, and the full 

solution using a computer program (10); radius ratios of 0.80 and 0,99 

“6 

are shown for values of the parameter, -K^/2 = 6x10 . The approximate 

formula (Eq, [39]) predicts a 27% decrease in mass leakage due to rota- 
tion for the externally pressurized mode of the compressor seal exampleo 
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Even though rotation has a substantial effect on mass leakage, the 
pressure distribution is affected very little* The opening force changes 
by less than 1/2% and there was virtually no change in the center of 
pressure* Hence in this case the seal pressure balance is not affected 
by rotation. 

The relative surface rotation appears to be negligible for many 
cases where the radial pressure differential is large, the gap is small, 
and the rotational speed is moderate. Such a case is evaluated in the 
appendix. This contrasts with the case for a liquid medium (incompres- 
sible fluid) , where the centrigual force in laminar radial flow can have 
a substantial effect because of the higher density. For example, water 
at room temperature is about three orders of magnitude more dense than 
air. The comparable incompressible viscous flow analysis (5) yields the 
following approximate formula. 


AM 

^^static 



[40] 


If a very small pressure differential exists across the sealing dam, 
the rotational effects on the radial flow are important even for moderate 
speeds. An example of such a case is also given in the appendix. 

The rotational flow is always important for determining the power 
loss (due to viscous shearing) and for determining the transition to 
turbulent flow. The rotational flow component of the velocity may be 
responsible for the flow becoming turbulent in the circumferential direc- 
tion (shear flow direction) , which will mean the entire flow field (radial 
pressure flow) will be turbulent c 
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Figure 7 illustrates a simplified envelope of possible sealing dam 
flow regimes 5 showing the region in which this analysis is valid* For 
the hydrostatic case^ the viscous flow analysis presented is valid for 
Rej^(h/AR) < 1* The parameter governing the flow regime for radial flow 
in most seal problems was found to be the Mach number rather than the 
pressure flow Reynolds number. This is true for very small gaps since 
calculations of (4) indicate that the flow can become choked in the 
radial direction before the transition (laminar to turbulent) pressure-- 
flow Reynolds number is reached. Choking of the flow, of course, will 
occur when M = 1. The limit for this analysis to be valid is M < l//y; 
beyond this, the isothermal assumption no longer holds. (This is shown 
in (4) . Reference (4) further shows that the adiabatic assumption is 
valid to M = 1.) 

For radial pressure flow only, transition to turbulent flow will 
occur by Tollmien-Schlichting waves. The minimum transition Reynolds 
number average velocity and hydraulic diameter) will 

be between 2300 and 6000; that is, between the minimum transition 
Reynolds number for fully developed pipe flow and what appears to be the 
theoretical minimum transition Reynolds number for plane Poiseuille 
flow (12) in infinitesimal disturbance hydrodynamic stability theory. 

The radial flow situation of the sealing dam approximates plane 
Poiseuille flow because the aspect ratio used in the hydraulic diameter 
calculation is the same. 

For the case of only circumferential shear flow with no imposed 
radial pressure gradient, the flow will remain laminar until a critical 
rotational Reynolds number is exceeded; then the flow will become turbu- 


20 



lent. For this narrow gap sealing dam analysis (AR >> h) , the critical 
rotational Reynolds number for transition appears to be simple Couette 
flow transition Reynolds number. Thus, Re^ = R^h/v = 1900 for tran- 
sition (13); the smaller the gap, the longer transition speed is de- 
layed. Thus, the wall has a stabilizing effect on the flow, and this 
should be considered in selecting the design gap (e.g,, for a 1-ft- 
diam seal operating at room temperature, the transition speed would be 
71 500 rpm for a 1-mil gap, but would be 179 000 rpm for a 0.4-mil gap). 

When there is both circumferential shear flow and radial pressure 
flow present, the net flow is a spiral flow. To the authors^ knowledge 
the transition from laminar to turbulent flow for this case has not 
been investigated. 

Most flow regime transitions have been experimentally observed for 
incompressible flow between a disk and a confined housing, but there is 
reason to believe that transition will also occur for a compressible 
fluid in this range of Reynolds number. Theodorsen and Regier (14) 
conducted an experimental investigation on a rotating disk in an infinite 
compressible fluid medium and found that laminar-turbulent transition and 
drag are independent of Mach number. The highest transition Reynolds 
number was 3x10^ for the most highly polished disk. However, 

experimental work is needed to find the transition Reynolds number for 

compressible flow between two disks separated by a very narrow gap in 

relative rotation which is representative of gas film seal operation. 

In this analysis, which is valid for a modified Reynolds number, 
Re^(h/AR) << 1, the circumferential velocity v can be found inde- 
pendently of the radial velocity u (i.e., v does not depend on u, but 
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u depends on v (see Eq« [13]). The radial velocity u is affected by 
the rotation through the centrifugal force component in the radial di- 
rection (see Eq. [12]). In a higher-order analysis, when Re^(h/AR) is 
no longer small, the fluid velocity component in the circumferential 
direction will depend on the radial pressure flow, which will, in turn, 
depend on the rotating flow velocity component. 

The region where the radial and the tangential flow are coupled in- 
volves the solution of two simultaneous nonlinear partial differential 
equations and is beyond the scope of this analysis. This more complex 
formulation would be necessary for a stability analysis and for a more 
detailed examination of the flow field. In turbulent flow, of course, 
both the circumferential and radial flow velocity components depend on 
each other. However, the simplified analysis presented herein is repre- 
sentative of many shaft face seal applications. 

In general leakage changes due to rotation must be considered with 
respect to other influences. Examples of other effects which under 
practical conditions may be of the same order or even more important 
are seal face deformations, surface wobbling, misalinement , and eccen- 
tricity. 

Rotational Effects on Pressure and Velocity Profiles 

The pressure distribution in the radial direction was found from 
Eq. [25] « The internally pressurized case is shown plotted in Fig. 8(a) 
for -K^/2 = 0.02. Also shown in Fig. 8(a) is the hydrostatic radial 
pressure distribution, which was found from Eq. [25] with Q - 0. Note 
the sealing dam force is slightly greater for rotation. The radial 
pressure distribution for the externally pressurized case is shown in 
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Fig. 8(b). Here the sealing dam force is slightly less than for the 
hydrostatic case. As previously discussed, these pressure distributions 
are only slightly different than the hydrostatic case. Hence for most 
practical purposes the pressure balance is unaffected. 

The radial velocity profiles for both internally and externally 
pressurized cases were found for ” 0.99, and 

-K^/2 = 0.02 and are shown in Figs. 9(a) and (b) . These radial velocity 
profiles were obtained from Eq. [16]. Also shown are the hydrostatic 
and centrifugal force velocity components. As expected, the centrifugal 
force component is effective in reducing the leakage in the externally 
pressurized case whereas the opposite occurs in the internally pres- 
surized case. The maximum velocity does not occur at the midgap position 
for the centrifugal force velocity component as it does for the hydro- 
static component. 

The shear flow (circumferential direction) velocity profile is the 
classical simple Couette flow profile. 

SUMMARY OF EESULTS 

The flow across a parallel sealing dam is analyzed for steady, 
laminar, subsonic, isothermal compressible flow with rotation. The 
assumed conditions are representative of many shaft face seal applica- 
tions. The analysis considered the effect of relative seal face rota- 
tion on mass leakage, pressure distribution, velocity distribution, and 
net pressure force; the hydrostatic case was used for comparison. The 
following pertinent results were obtained: 

1. Rotation has significant effects on seal leakage flow at high 
speeds. (The seal example analyzed has a 42% decrease in calculated 
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leakage for an externally pressurized sealing mode*) The seal pressure 
balances however^ appears to be substantially unaffected by rotation* 

2* When the radial pressure differential is very small j rotational 
effects are also important at moderate speeds* (The seal example ana- 
lyzed has a 17-percent increase in calculated leakage for a 0*2-psi 
pressure differential at a moderate speed*) 

3* A hydrostatic analysis of the radial pressure flow is a good 
approximation of the rotating seal under the condition that the ratio 
of the reference radial pressure flow velocity to reference rotational 
velocity is greater than the square root of the modified Reynolds 
number . 


4. The error in estimating the mass leakage by using the hydrostatic 
radial flow formula was analyzed, and a pertinent error analysis equation 
is given, 

5. The region of validity of the analysis was defined as follows; 

a* The analysis is valid until transition of the circumferen- 
tial shear flow to turbulence. 

b. The analysis is valid until the radial pressure flow 
approaches a Mach number of l/Zy^ where y is the ratio of specific 
heats, or until the radial pressure flow becomes turbulent, 

6. The circumferential shear flow is important in that at high 
rotative speeds it can cause the flow field to become turbulent. 
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APPENDIX - EXAMINATION OF CENTRIFUGAL FORCE EFFECT 
(ROTATION) ON RADIAL AIRFLOW 

Consider the following case. Let v = 2x10 ^ ft^/sec, T = 500° R, 

R = 111 foot, ^ - 2x10° ft^/sec - R, N = 2300 rpm (this can be consid- 
ered a moderate speed), and h = 0.001 inch. Hence, Re^ = 25, much less 
than the critical transition Reynolds number for simple Couette flow, 
which is about 1900. Since |k^/ 2| = 4x10 > Eq. [39] can be used 

to find the change in mass leakage caused by rotation. 

Case I - Very Small Pressure Differential 
When 

P^ = 100.2 psia, R^ = 5.70 in.^ _ 

> at R = 1/2 ft 

= 100 psia, R2 = 6.30 in, J 

Eq. [39] yields ^/^static “ percent. 

Hence, rotation also has a substantial effect on mass leakage when 
a very small pressure differential exists. In this example, the radial 
flow Mach number is very small, and the incompressible formula (Eq. [40]) 
predicts about a 25-percent change. 

Case II - Moderate Pressure Differential 
When 

P^ = 100 psia, R^ = 5,70 in. 

P 2 = 150 psia, R 2 = 6.30 in. 

Eq. [39] yields 
AM . 

T s -0.1 percent 

^static 

Thus, the above example illustrates that centrifugal force usually does 
not have an appreciable effect on air mass leakage for moderate and large 
pressure ratios and moderate rotational speeds* 
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Figure 3. - Model of sealing dam. 
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(A) SEALING DAM MODEL USED IN ANALYSIS. MODEL CHARACTERISTICS: 
NO AXIAL FLOW; NARROW GAP, h«AR. 

(B) RADIAL DIFFUSER MODEL. MODEL DOES NOT CONFORM TO MODEL OF 
FIGURE 4(A) WI€N LARGE AXIAL aOW RATES EXIST AND/OR AR~ ^(h). 

(Cl HYDROSTATIC OR EXTERNALLY PRESSURIZED THRUST BEARING MODEL 
MODa DOES NOT CONFORM TO MODEL OF FIGURE 4(A) WHEN SIGNIFI- 
CANT AXIAL FLOW THROUGH AN ORIFICE EXISTS. 

(D) RADIAL STEP SEAL MODEL. MODEL OF FIGURE 4(A) IS VALID WHEN 
hi» h, (THEN Pi2 = Pi) AND INNER CAVITY VaOCITY IS NOT 
SIGNIFICANT. 
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Figure 5. - Effect of rotation on mass leakage 
for a proposed space power compressor seal 
(11) for both internally and externally 
pressurized modes of operation. 


Figure 4. - Various seal and bearing models illustrating conditions of analysis 
validity. 
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Figure 6. - Change of mass flow due to rotation as function 
of pressure ratio. Comparison of numerical solution with 
approximate formula. 
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Figure 7. - Envelope of possible sealing dam flow regimes, 
illustrating region of analysis validity. 
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(A) INTERNALLY PRESSURIZED 
CASE. PRESSURE RATIO, 
Pl/P2= 1.2. 


(B) EXTERNALLY PRESSURIZED 
CASE. PRESSURE RATIO, 
Pj/P2 = 0.833. 

Figure 9. - Radial velocity pro- 
files. moving wall in 9-plane. 
Radius ratio, Rj/R 2 “ 0.99. 
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Figure 8 , - Effect of rotation on 
radial pressure distribution. 



